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Magnetic reconnection occurring in collisionlcss environments is a multi-scale process involving both ion and 
electron kinetic processes. Because of their small mass, the electron scales are difficult to resolve in numerical 
and satellite data, it is therefore critical to know whether the overall evolution of the reconnection process 
is influenced by the kinetic nature of the electrons, or is unchanged when assuming a simpler, fluid, electron 
model. This paper investigate this issue in the general context of an asymmetric current sheet, where both 
the magnetic field amplitude and the density vary through the discontinuity. A comparison is made between 
fully kinetic and hybrid kinetic simulations of magnetic reconnection in coplanar and guide field systems. The 
models share the initial condition but differ in their electron modeling. It is found that the overall evolution 
of the system, including the reconnection rate, is very similar between both models. The best agreement 
is found in the guide field system, which confines particle better than the coplanar one, where the locality 
of the moments is violated by the electron bounce motion. It is also shown that, contrary to the common 
understanding, reconnection is much faster in the guide field system than in the coplanar one. Both models 
show this tendency, indicating that the phenomenon is driven by ion kinetic effects and not electron ones. 



I. INTRODUCTION 

Magnetic reconnection is a universal plasma phe- 
nomenon enabling large scale plasma structures frozen 
in the magnetic field to change their magnetic connec- 
tivity while transferring a substantial part of the mag- 
netic energy stored in current sheets into the thermal 
and bulk kinetic energy of the surrounding plasma 1 - 7 . In 
most systems where reconnection is thought to play a key 
role, the amount of magnetic flux reconnected per time 
unit, the so-called reconnection rate, appears as a critical 
parameter for their large scale evolution. Understand- 
ing which physical mechanism controls it, and at which 
scale, is therefore an important issue from the funda- 
mental and modeling viewpoint. In collisionless systems, 
the vastly different inertia of electrons and ions makes 
magnetic reconnection a multiscale process^ 7 . Because 
it creates accelerated and heated flows and macroscopi- 
cally changes the transport of the plasma in magnetized 
systems, magnetic reconnection can be thought as a fluid 
phenomenon. However, the collisionlcss nature of the 
systems in which it occurs involves kinetic processes for 
which the importance regarding this large scale evolu- 
tion is still poorly understood. Although the question 
really concerns both species, knowing to what extent the 
kinetic behavior of the electrons plays a role in control- 
ling the overall dynamics is particularly critical to un- 
derstand, observe/measure and model the reconnection 
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process, since it concerns scales that are, from all those 
viewpoints, very difficult to resolve, in comparison to ion 
scales. 

Over the years, several studies have addressed this is- 
sue and different conclusions have been reached. On one 
hand, electron kinetic physics is seen as the fundamen- 
tal way to break the flux freezing constraint in collision- 
less systems, but is thought not to affect the global pro- 
cess, which is rather controlled by the coupling of ion 
dynamics and electron fluid physics 3 -* 4 -. On the other 
hand, electron kinetic physics is seen as a key ingredient, 
which can make the whole process unsteady by enabling 
the frequent production of magnetic flux ropes£~— . The 
applicability of our understanding is, however, most of 
the time limited by the assumption of an initial symme- 
try across the current sheet, which might be adequate 
for the modeling of some environments like the Earth 
magnctotail, but seems oversimplified in most other cur- 
rent layer systems like the Earth magnctopause, the so- 
lar wind and solar coronal loops. The few numerical 
studies which focused on kinetic reconnection in asym- 
metric current sheet, have revealed a great number of 
new features challenging our current understanding of 
magnetic reconnection^—. Among them, recent stud- 
ies have revealed a rich variety of electron scale kinetic 
processe a 10 ' 11 ' 13 ' 14 ' 16 . Surprisingly, many of them were 
not confined to the reconnection site as one could have 
expected, but rather expand over large regions of the re- 
connection exhaust. Their omnipresence raises the ques- 
tion of whether including electron kinetic physics is a 
mandatory requirement for an adequate modeling of large 
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scale reconnection of magnetic field in these more gen- 
eral systems, and to what extent neglecting it changes 
the overall dynamics. Along the same lines, the initial 
condition of kinetic simulations of asymmetric reconnec- 
tion models is never at a kinetic steady state but in a 
fluid equilibrium, which very rapidly evolves toward a 
self-consistent quasi-steady state with an internal struc- 
ture very different from the prescribed one^i. The role of 
kinetic electrons in the establishment of such structure is 
unknown and it is important to understand to what ex- 
tent neglecting this physics impacts the structure of the 
current sheet where reconnection develops. 

In this paper, we focus on the role of electron kinetic 
physics in asymmetric magnetic reconnection. To do so, 
we present the results of two kinetic simulations sharing 
the same initial configuration but differing in the model- 
ing of the electron population. One code solves the com- 
plete Vlasov-Maxwell system, and is commonly referred 
as a fully kinetic Particle-In-Cell (PIC) codeii, while the 
other one assumes the electrons behave as an isother- 
mal fluid without bulk inertia and the ions are treated as 
particles, and is usually referred as a hybrid PIC codei^. 
These latter assumptions are of course not generally sat- 
isfied in real systems, but their consequences regarding 
the reconnection process are not well understood. The 
second part of this paper is dedicated to the description 
of the simulation methods with their respective physical 
assumptions. The third part describes the initial condi- 
tions used in this work and investigates to what extent 
the kinetic nature of electrons impacts the initial internal 
structure of the current sheet after the re-configuration of 
the initial fluid equilibrium. The fourth part presents the 
reconnection rate obtained from both simulations. The 
fifth part discusses the structure of the current sheet and 
the sixth and last section summarizes and discuss our 
findings. 



II. NUMERICAL MODELS 

This section describes the numerical models used in 
this study. The fully kinetic model treats both electrons 
and ions as a collection of particles feeling the electric E 
and magnetic B fields. The force they feel is described 
by eq. (JXJ) , where m s , q s and v s are the mass, charge and 
velocity of the particle of species s. The electromagnetic 
fields are obtained by solving the Maxwell-Ampere © 
and Maxwell-Faraday (|3|) equations, where j is the total 
electrical current density, c is the speed of light in vac- 
uum and fio the magnetic permeability of vacuum. The 
details of the algorithms have already been explained in 
a previous paper—. For the runs presented in this paper, 
we have used an electron to ion mass ratio m e /m; = 1 /25 
and the ratio between the plasma frequency and the elec- 
tron cyclotron frequency is uj pe /uj ce = 4. 
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= V x E 



(1) 
(2) 
(3) 



In the hybrid code, only the ions are treated as a col- 
lection of particles. Electrons, on the other hand, behave 
as a fluid, i.e. no kinetic behavior is included in their 
modeling. Furthermore, because of their small mass com- 
pared to ions, we also neglect quasi-neutrality scales and 
assume the electron density equals the ion one. Consis- 
tently, the displacement current is neglected in Maxwell- 
Ampere equation ([2]). The electron bulk inertia is also 
neglected. Their fluid momentum equation is thus used 
to update the electric field, based on these approxima- 
tions (dJ), where v e is the electron fluid velocity, P e their 
scalar pressure and R a dissipation term. 
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x B VP e 
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(4) 



In this study, we choose an hyperrcsistivity to model 
the dissipation, R = — i/V 2 j, where v = 5 10~ 4 . 
Such dissipation is analogous to an electron viscosity 
mechanism, which is coincidcntally similar to what 
occurs at the reconnection site in symmetric fully 
kinetic models^. Its implementation in the modeling 
of asymmetric magnetic reconnection has furthermore 
been recently shown to have important consequences 
regarding the overall evolution of the system^!. Details 
of the algorithms used in the hybrid code can be found 
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Both models are have 2.5 dimensions. The plane in the 
fully kinetic model is called x — z whereas in the hybrid 
model it is called x—y. The domain is considered periodic 
in the x direction and closed by perfect conducting walls 
in the other one. The data presented in this paper is in 
normalized units. The magnetic field is normalized to an 
arbitrary magnetic field Bq, the density to an arbitrary 
density uq, the time to the ion cyclotron period based 
on the ion mass and charge and on the magnetic 
field Bq. The lengths are normalized to the ion inertial 
length Si based on the density no. The electric field is 
therefore normalized to VaBq, where Va = -Bo/V n o TO iMo 
is the characteristic Alfven speed in the system. 



III. INITIAL CONDITION AND EQUILIBRIUM 

In this section, we describe the way the system stud- 
ied in this paper is initialized. To simplify the problem, 
we assume the current sheet to be a one dimensional 
tangential discontinuity where the plasma properties and 
electromagnetic fields vary between two asymptotic and 
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BJB a 0.5- 




FIG. 1. The solid line represents the initial in-plane magnetic 
field profile given by eq. 0. The dashed line represents the 
initial density given by eq. ©. 



uniform values. We also assume the system to be in 
steady state before reconnection occurs. Following these 
hypotheses, the density n, temperature T and magnetic 
field B profile must satisfy the ID pressure balance con- 
dition ([5]), where ks is the Boltzmann constant. 



nksT H = est 

2/io 



(5) 



We choose the following profiles for the density ^ , the 
magnetic field ([7]) where A = 0.5. The ion temperature 
Ti is obtained from equation ([5]) and T = Tj +T e , where 
T e = 0.2T; is the electron temperature. 



i(y) = l-l(s(y) + S(yf) 



B x (y) = - + S(y) 



S (y) = tanh 



y - ya 
A 



(6) 
(7) 
(8) 



This initial condition has been used in previous 
studie a n ' 13 i 24 and is illustrated on Fig. [TJ For the sake 
of generality, we will consider two canonical configura- 
tions : the first includes an out-of-plane component of 
the magnetic field, which, for simplifying purpose, is as- 
sumed to be uniform and equal to B g f = 1 and the sec- 
ond one is a coplanar current layer (B g f = 0). Equa- 
tions ([6][7]) arc defined in the hybrid coordinate system, 
replacing y by z defines them in the plane of the fully 
kinetic model. The domain is 64<5i in the downstream di- 
rection and 25.6<5i in the upstream direction, j/o denotes 
half of the upstream size. Because our system is peri- 
odic in the x direction, the tearing instability developing 
within the current layer would eventually saturate with a 



fixed number of magnetic islands instead of ejecting them 
away and form a dominant X line where plasma jets are 
formed. To prevent this artificial evolution, and because 
we are mostly interested in the non-linear quasi-steady 
reconnection regime, we initially trigger magnetic recon- 
nection with a local magnetic perturbation of amplitude 
SB, centered in the middle of the domain. To isolate 
the effect of the initial perturbation on the system from 
the one originating from lack of electron kinetic physics, 
the hybrid calculations are repeated with a perturbation 
twice as large as the one used in the original runs. The 
simulations presented in this paper are summarized in 
the table HI 

TABLE I. Summary of the simulations presented in this paper 



Run 


Model 


B 9f 


SB 


FG 


full PIC 


1.0 


0.1 


FC 


full PIC 


0.0 


0.1 


HGi 


Hybrid PIC 


1.0 


0.1 


HCi 


Hybrid PIC 


0.0 


0.1 


HG 2 


Hybrid PIC 


1.0 


0.2 


HC 2 


Hybrid PIC 


0.0 


0.2 



The particles are initially loaded in Maxwellian dis- 
tribution functions which locally have the density and 
temperature described above. Although it satisfies the 
pressure balance condition ([5]), this initial condition is 
not a solution of the steady Vlasov-Maxwell system. As a 
result finite Larmor radius effect will modify the internal 
structure of the current laye r 11 ! 25 where both electrons 
and ions find a new and self-consistent force balance. 

From the modeling and fundamental viewpoint it is 
important to understand to what extent the kinetic na- 
ture of electrons participates to the establishment of this 
new equilibrium^ - — , and what are the consequences of 
ignoring this part of the physics 2 ^. A good test is to 
compare the early evolution of the system between fully 
kinetic and hybrid kinetic simulations. Figures [5] and 
[3] show the ion and electron force balance at t = 5 in 
both codes, in a cut in the cross current direction, for 
the coplanar and guide field configurations, respectively. 
The curves are obtained after an average of the force 
profiles between x = and x = 15 to reduce the noise. 
As one can notice, the amplitude and variations of the 
forces are very similar in the fully kinetic and hybrid 
kinetic systems, for ions as well as for electrons, and for 
the coplanar and guide field cases. In the coplanar case, 
the ion equilibrium is mostly the result of a balance 
between the pressure force — VP; and the electric force 
enE, while in the guide field case the magnetic force 
n p9p v i x B becomes as important as the electric force 
in balancing the pressure gradient force, indicating that 
ions are more magnetized. The electron equilibrium 
however, mostly concerns the balance between the 
electric force n e <7 e E and the magnetic force n e q e v e x B, 
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and that in both the coplanar and guide held case. 
In a movie, available online, one can notice that the 
time scale over which the layer oscillates and emits 
waves are also identical in the hybrid and fully kinetic 
runs. As a result of this remarkable ressemblance, one 
can conclude that the self-consistent fluid equilibrium 
found by the system is mostly a consequence of the 
kinetic behavior of ions and not of the electrons, for 
which one can safely neglect the kinetic nature. Systems 
having an initially non-uniform temperature might show 
more differences and this have to be checked in future 
models. Let us notice, however, that the current sheet 
has initially a sub ion scale thickness, which is quite 
small compared to the Earth magnctopause current 
sheet for instance, and that the mass ratio used in this 
study is unrcalistically large. As a consequence of these 
two effects combined, our model largely emphasizes the 
role of electron finite Larmor radius effects, if any. In 
a realistic configuration, considering the electrons as 
a fluid should therefore be an even better approximation. 
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FIG. 2. Force balance across the current sheet for the ions 
(top panels) and the electrons (bottom panels) at t = 5 in the 
fully kinetic (left panels) and hybrid kinetic (right panels) 
simulations initialized with the coplanar configuration. 
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FIG. 3. Force balance across the current sheet for the ions 
(top panels) and the electrons (bottom panels) at t — 5 in the 
fully kinetic (left panels) and hybrid kinetic (right panels) 
simulations initialized with the guide field configuration. 



IV. RECONNECTION RATE 

In this section, we compare the reconnection rate ob- 
tained from the hybrid and fully kinetic models. The 
reconnection rate is the amount of magnetic flux recon- 
nected per time unit. In a two-dimensional configuration, 
this definition corresponds to the electric field compo- 
nent perpendicular to the reconnection plane, measured 
at the X line. Technically, this electric field can be ei- 
ther directly measured in the simulation or be calculated 
as the time derivative of the out-of-plane component of 
the vector potential <f> at the X line, which in 2D repre- 
sents the variation of the reconnected flux. The latter 
method is preferred since it involves much less noise in 
the results than the former. Figure 2] shows the recon- 
nection rate for both numerical models in the coplanar 
and guide field configurations. The X line itself is found 
as being the saddle point of the out-of-plane component 
of the vector potential and the measure is averaged over 
an area of 0.04<5^. The top panels show the time evolu- 
tion of the reconnection electric field. The bottom panels 
show this rate as a function of the reconnected flux. All 
plots share the same tendency, where the rate first in- 
creases over a certain time interval, reaches a maximum 
value and slowly begins to decrease. This late phase is 
the consequence of the finite size of the domain, which 
limits the region over which the reconnected flux can be 
expelled and the amount of magnetic flux that can be 
reconnected. 

Let us first look at the guide field configuration. On 
the top right panel, the hybrid and full PIC simulations 
sharing the same perturbation have a substantially differ- 
ent time evolution of the reconnection rate. However, it 
is worth noticing that their maximum reconnection rate 
is quite similar, the main difference being in the time 
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FIG. 4. Top panels : Reconnection rate as a function of time for the coplanar (left) and guide field configurations (right). 
Bottom panels : Reconnection rate as a function of the reconnected flux <f) for the coplanar (left) and guide field configurations 
(right). On all panels, the blue (dash-dot) and red (solid) curves are respectively obtained from the fully kinetic and hybrid 
runs sharing the same perturbation amplitude. The green curve (dash) is obtained from a hybrid run having a perturbation 
with an amplitude twice larger. 



period required to reach it. For the same perturbation 
amplitude, the hybrid model requires more time to reach 
the maximum rate than the fully kinetic model. Doubling 
the perturbation amplitude makes both models almost in 
phase. Because the two models differ only by the han- 
dling of the electron physics, these results suggest that 
the electron kinetic physics in one case, and the hypcr- 
rcsistivc dissipation in the other case, do not respond the 
same way to the initial perturbation. This comes from 
the fact that the early evolution of the process consists 
in a local collapse of the magnetic reversal down to the 
dissipation scale where convection and dissipation are in 
competition. Because the dissipation mechanism is dif- 
ferent in both models, the time required to build the 
current sheet is different. In a fully kinetic simulation, 
one might expect this phase to depend on the electron 
to ion mass ratio. Fully kinetic codes are also subject to 



a larger noise than hybrid codes, which can also accel- 
erate the development phase of the process. Looking at 
the top panel for the same configuration reveals that the 
reconnection rate as a function of the reconnected flux is 
identical for all simulations. Consequently, the role of the 
perturbation is just a time lag on the reconnection rate, 
which, fundamentally, rather depends on the phase of the 
process, i.e. on the upstream properties of the field and 
plasma, regardless of the kinetic nature of the electrons. 

The situation is a bit different for the simulation of 
the coplanar configuration. Here again the rate of the 
fully kinetic model increases faster than the one of the 
hybrid model, however the time lag of the latter is not 
as much as in the guide field case. As a result, doubling 
the perturbation amplitude, which, again, shortens the 
development phase, makes the hybrid model to reach the 
maximum before the fully kinetic one. As for the guide 
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field scenario, the electron/dissipation physics does not 
respond the same way to the perturbation. However, con- 
trary to the guide field case, it is unlikely that one finds 
a perturbation amplitude for which both models would 
be in phase with the same amplitude. Indeed, while they 
do not differ drastically, the maximum rate of the fully 
kinetic model is significantly larger than the hybrid one, 
which, itself, appears moreover not to depend on the am- 
plitude of the perturbation. When looked as a function 
of the reconnected flux, the rate of both hybrid simula- 
tions appears to be more similar and separated from the 
one of the fully kinetic models. These observations sug- 
gest that in the coplanar case, the reconnection rate is 
sensitive to the nature of the electron physics, and that 
kinetic electrons apparently speeds up the process. 

Finally, it is worth noticing the difference between 
the coplanar and guide field configurations. One can 
clearly see that the guide field runs have much larger 
reconnection rates than the coplanar ones. While guide 
field reconnection is as fast or slower than antiparallcl 
reconnection in symmetric systems because in the lat- 
ter the plasma is more compressible, our results suggest 
that this well accepted scenario cannot be extrapolated 
to asymmetric systems for which a better explanation 
has to be found. Because asymmetric current sheets are 
much more prevalent than symmetric ones, this result 
has broad consequences, in particular if one considers 
magnctopause reconnection. This finding will be stud- 
ied in detail in a forthcoming paper. Let us remark at 
this point that the effect is seen in both hybrid and fully 
kinetic models, indicating that ions are the main contrib- 
utor to this effect and not electrons. 



V. STRUCTURE OF THE CURRENT SHEET 

During the reconnection process, the magnetic field 
reversal at the X line itself is associated with a current 
sheet mainly supported by the lightest species in the 
plasma, namely the electrons. At this location, the 
current density is sustained through the acceleration 
of the electrons by the reconnection electric field and 
dissipated as the hot and fast electrons leave the layer 
by recoupling to reconnected field lines. A steady 
state is reached when these two effects balance each 
other—. From the fluid viewpoint, this region of space 
is characterized by the competition between advection 
and dissipation mechanisms, and it is not clear to what 
extent the underlying kinetic behavior of the electrons 
controls the structure of the current sheet. Figure [5] 
shows the out-of-plane current density for the hybrid 
and fully kinetic simulations in the coplanar and guide 
field configurations at t = 35. Figure [5] represents a slice 
of the current density, the ion current density and the 
electron current density at the same time for the same 
simulations, along the upstream direction and through 
the X line. Both figures show the same features overall: 
the current sheet has a sub-ion scale mainly supported 



by an electron current, as expected. Looking at Fig. [5] 
and in particular at the right panels, we can see that 
despite the difference in the symmetry originating from 
the different coordinate systems and the same guide 
field sign, both simulations look very similar from the 
ion scale down to the current sheet scale. Both models 
show a left right asymmetry that is a consequence of 
the initial guide field together with the Hall effect. The 
negative current density is strongly enhanced at the 
reconnection site and spreads on the top separatrices 
with more pronounced values on the big island side 
of the X line. Both models also show a current layer 
with an opposed sign on the bottom separatrix with 
values more pronounced on the smaller island side. The 
amplitude of the current density at the X line itself is 
roughly similar, but this will probably depend on the 
mass ratio in the fully kinetic model and on the value 
of the hyperresistivity in the hybrid kinetic model. The 
bottom panels of Fig. [5] reveal the very similar structure 
of the total current density between the two models. 
In both cases, the current sheet docs not consists of 
a single peak but has a smaller peak of opposite sign 
on the weak field side. In both cases, the respective 
part of the current density supported by the ions and 
the electrons looks similar in amplitude and the spatial 
extent of the non-zero current density is also identical. 
One can notice that gradients look slightly sharper in 
the hybrid results than in the fully kinetic ones, this 
effect is probably controlled by the electron to ion mass 
ratio and the value of the hyperresistivity. 

As for the reconnection rate, the coplanar runs are 
again revealing more differences between the hybrid and 
fully kinetic models. The left panels of Fig. [5] show that 
the current sheet in the hybrid run is much more local- 
ized than its fully kinetic counterpart. In the fully kinetic 
simulation, the current sheet is both broader and longer 
in the upstream and downstream directions, respectively. 
As a result, the current density is much weaker than in 
the hybrid run. Looking at the top panels of Fig. [51 
one can see that the ion current density is very similar 
between both runs. It has the same amplitude and its 
variations are identical. The total current density is 
however different. While the hybrid one has roughly the 
same structure as in the guide field configuration, with 
the peak of opposite sign on the weak field side, the fully 
kinetic current density looks more symmetric, even if 
one can notice a slight increase on the weak field side too. 

The fact that the structure of the hybrid current sheet 
differs from the fully kinetic one mainly in the coplanar 
case is a consequence of the lack of electron confinement 
in the vicinity of the X line. Like the ions, the electrons 
see a vastly different field on both sides of the current 
sheet, and they are less confined on the weak field side 
than on the strong field side. The guide field amplitude 
is large enough to change the ion magnetization in the 
current sheet, therefore it is also sufficient to magnetize 
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FIG. 5. Out-of-plane current density at t = 35 for the fully kinetic (top panels) and hybrid runs (HG2 and HCi, bottom 
panels) in the coplanar (left panels) and guide field (right panels) configurations. The current density of the fully kinetic model 
as been multiplied by —1 so that it can be compared more easily with the hybrid results obtained in the different coordinate 
system. 



the electrons, which are much lighter. As a result, the 
guide field scenario is less sensitive than the coplanar 
one to electron kinetic effects resulting from non-local 
mixing of populations. In the hybrid coplanar simulation 
however, electrons physics is local, which therefore leads 
to a more confined current sheet. Due to the unrealistic 
electron mass used in the simulations, we expect that 
lighter electrons will result in their better confinement 
and less differences between the hybrid and fully kinetic 
models. 



VI. MOTION OF THE X LINE AND DISSIPATION 
REGION 

Our hybrid model does not include non-gyrotropic, or 
even anisotropic electron pressure. It is therefore inter- 
esting to tests to what extent, the features for which these 



effects have been identified as a key ingredients, differ 
between the hybrid and fully kinetic runs. We choose to 
investigate two of those features. The first one is the ex- 
tent of the dissipation region around the X line. It has re- 
cently been shown^S that the electron scale current layer 
was not necessarily a good proxy of a dissipative process. 
Another recent study^l then proposed a scalar quantity 
D e to measure non-ideal energy transfers from the elec- 
tromagnetic fields to the plasma. In a collisionless en- 
vironment, irreversible energy transfers are the result of 
complicated mixing is phase space that macroscopically 
appear in the non-gyrotropic components of the electron 
pressure tensor. In the hybrid model, these terms are 
missing, however, they are modeled by a simple uniform 
hyperresistivity. Fig. [7] shows the dissipation measure 
D e in the hybrid and fully kinetic runs, for the coplanar 
and guide field configurations. D e is strongly localized in 
the reconnection region in both hybrid and fully kinetic 
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FIG. 6. Out-of-plane current density through the X point at 
t — 35 for the fully kinetic and hybrid runs (HG2 and HCi) in 
the coplanar and guide field configurations. For all panels, the 
total current density is represented in black, the ion current 
density in red and the electron one in blue. As in Fig [5] the 
current densities in the results obtained from the fully kinetic 
code have been multiplied by —1 to ease the comparison. 



runs, and is better confined in the hybrid runs, since no 
electron finite Larmor radius effect is occurring. We can 
again notice the confinement effect the guide field has 
on electrons, as the dissipation region is, in this config- 
uration, better localized than in the coplanar case. One 
puzzling result is that the strongest dissipation seems to 
be shifted from the actual position of the X line, it ap- 
pears to be collocated with the maximum of the current 
density (Fig. [5]). This could have already been notice 
in a previous worked were the same initial condition was 
used. These features are seen in fully kinetic and hybrid 
runs, which indicates that they are not simply resulting 
from an electron kinetic effect. A detailed investigation 
of the structure of the dissipation region will be the topic 
of a future study. 

The second feature we will discuss is the motion of 
the X line often seen in asymmetric configuration with a 
guide field and often associated to an electron diamag- 
netic effect^, i.e. an electron pressure effet. As one can 
see in the figures [5] and [7j the X line in the fully kinetic 
model is approximately located at x = 33.2, whereas 
the hybrid X line appears closer to its original location 
(x = 32). To investigate whether the two models have 
a qualitative difference regarding the motion of the X 
line, we show, on FigjSl the position, in the reconnection 
plane, of the X line, in the hybrid and fully kinetic mod- 
els. Noticing that the apparent symmetry is a result of 
the different coordinate system, one can see a very strong 
similarity in the motion of the X line. First, the X line 
does not move very far from its initial location. Then, 
and more surprisingly, it begins to move in one direction 
and then starts moving in the other one until the end 
of the calculation. The effect is quantitatively seen in 



both hybrid and fully kinetic models. Small quantitative 
differences can be seen, the most important being that 
the fully kinetic X line travels farther than the hybrid 
one before changing direction a little bit later. There is, 
consequently, a slight timing difference for the reversal 
time (t w 27 in the fully kinetic run and t sa 20 in the 
hybrid one), which results, for the hybrid run, in an X 
line being coincidentally close to its original position at 
the time (t = 35) when Fig. [5] and [7] arc made. In both 
models, the X line also moves upward, with a slightly 
faster speed in the hybrid model. It is not clear why 
the X line changes its direction. The electron and ion 
diamagnetic drift speeds (not shown) are in opposite di- 
rections all along the simulation. However, their values 
around the X line exceed by far the velocity of the X line 
motion. It is possible that the first phase corresponds to 
an influence of the initial perturbation while the second 
phase is more a self-consistent feature of the reconnec- 
tion process. A more detailled investigation is beyond 
the topic of the present paper, which aims to focus on 
the overall evolution of asymmetric reconnection. 



VII. SUMMARY AND DISCUSSION 

In this paper we have compared hybrid and fully 
kinetic simulations of asymmetric magnetic reconnection 
to investigate the role of the kinetic nature of the 
electrons in the overall evolution of the system. We have 
chosen two configurations, one with a coplanar current 
sheet and a second with a uniform guide field. The 
initial condition is a fluid pressure balance with locally 
Maxwellian distribution functions. This initial state 
being not a Vlasov equilibrium, finite Larmor radius ef- 
fects rapidly change the internal structure of the current 
sheet and establish a new self consistent force balance 
while waves propagate away. We have shown that the 
force balance found by the system is identical in the 
full PIC and hybrid models, in amplitude, spatial and 
temporal scales, indicating that the process is controlled 
by the kinetic behavior of the ions and not by kinetic 
electrons. This should also be inspected and confirmed 
in other systems, having an initial temperature gradient 
for instance, however the very small thickness of our 
initial current sheet and the large electron mass used 
in this system indicate that the kinetic nature of the 
electrons might not be critical in real systems like the 
quiet magnctopause where gradients are usually at the 
ion scale. 

We have also shown that the reconnection rate 
obtained from the hybrid and fully kinetic results are 
fairly similar. One important difference is the response 
of the two models to the initial perturbation, the 
hybrid model taking more time to reach the maximum 
reconnection rate than the fully kinetic model. This 
effect originates from the distinct dissipation mechanism, 
which plays a critical role in the establishment of the 
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FIG. 7. Dissipation measure^ D e calculated at t = 35 for the fully kinetic (top panels) and hybrid (bottom panels) in the 
coplanar(left panels) and guide field (right panels) configurations. Notice the color range is adjusted for each panel. The small 
blue circle denotes the position of the X point, localized as the saddle point of the magnetic flux function. 



current sheet at the reconnection site and also possibly 
of the electrostatic noise of fully kinetic codes larger 
than in hybrid codes. Overall, the initial perturbation 
has been shown to affect the time dependance of the 
reconnection rate but its dependance on the phase of 
the process, i.e. on the upcoming magnetic field and 
plasma properties, stays unchanged. When plotted as 
a function of the reconnected flux, the difference in 
the reconnection rates between hybrid and fully kinetic 
models for the guide field configuration is negligible. 
However, the coplanar runs show some differences: the 
rate in the fully kinetic model being somehow larger 
than the hybrid one. These similarities and differences 
of the hybrid and fully kinetic models for the guide 
field and coplanar configurations, respectively, have also 
been shown in the structure of the out-of-plane current 
density. If it is very similar for the guide field case, it 
is appreciably different in the coplanar case. The guide 
field configuration leads to more confinement of both 
species inside the current layer, in particular of electrons. 
However this confinement is a kinetic process due to the 



mixing and the bouncing of particles inside the magnetic 
reversal, it can therefore occur for ions in both models 
but for electrons only in the fully kinetic model. As a 
result, without a guide field the ion current density gets 
broader in both models but the electron current density 
stays localized in the hybrid one, which changes the 
total current sheet structure. We expect the difference 
between the two models to diminish as the electron mass 
gets smaller, which should be tested in future studies. A 
last important point is the difference of the reconnection 
rate between the guide field system and the coplanar 
one. The guide field case is considerably faster than 
the coplanar system, which is surprisingly opposed to 
the common understanding of the effect of a guide field 
on reconnection. The fact that the hybrid and fully 
kinetic models both show the same tendency, to the 
minor differences explained above, indicates that ions 
are the primary responsible for this effect. This finding 
has broad consequences, for instance in the debate of 
whether magnctopausc reconnection would prefer guide 
field or coplanar configurations, and will be analyzed in 
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FIG. 8. Position Ax of the X line in the reconnection plane 
for the hybrid and fully kinetic model with respect to its ini- 
tial position. To ease the comparison, the blue curve shows 
the mirror with respect to x — from the actual position 
obtained from the hybrid model, which, otherwise be the op- 
posite because of the different coordinate system used. 



detail in a forthcoming paper. 

As a last step, we have studied the extent of dissipa- 
tive energy transfers and the motion of the X line, both 
effects being, associated to the electron pressure tensor. 
We have shown that the dissipation measure proposed 
recently^ is well localized in the reconnection region in 
both hybrid and fully kinetic models. In all models, the 
structure is slightly shifted from the actual X line and ap- 
pears collocated with the maximum of the current den- 
sity. In the presence of a guide field, the X line has 
been observed to move slowly although its motion has 
not been found related to an ion or electron diamagnetic 
effect. These results emphasize the need for a detailed 
investigation of the different processes occurring in the 
vicinity of the X line, which appear to differ substantially 
from the present understanding inferred from symmetric 
models. Future studies should also consider larger sys- 
tems and longer simulations times as they might possibly 
reveal unexpected behaviors as it did for reconnection in 
symmetric systems. 
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